On the Power of One Bit of Quantum Information 



E. Knill 1 , R. Laflamme 2 
1 MS B265, 2 MS B288, Los Alamos National Laboratory, Los Alamos, NM 87455 

(January 1998) 



In standard quantum computation, the initial state is pure 
and the answer is determined by making a measurement of 
some of the bits in the computational basis. What can be 
accomplished if the initial state is a highly mixed state and 
the answer is determined by measuring the expectation of a z 
on the first bit with bounded sensitivity? This is the situ- 
ation in high temperature ensemble quantum computation. 
We show that in this model it is possible to perform interest- 
ing physics simulations which have no known efficient classical 
algorithms, even though the model is less powerful then stan- 
dard quantum computing in the presence of oracles. 

Recent discoveries show that quantum computers can 
solve problems of practical interest much faster than 
known algorithms for classical computers . This has 
lead to widespread recognition of the potential benefits of 
quantum computation. Where does the apparent power 
of quantum computers come from? This power is fre- 
quently attributed to "quantum parallelism" , interfer- 
ence phenomena derived from the superposition princi- 
ple, and the ability to prepare and control pure states ac- 
cording to the Schrodinger equation. Real quantum com- 
puters are rarely in pure states and interact with their 
environments, which leads to non-unitary evolution. Fur- 
thermore, recent proposals and experiments using NMR 
at high temperature to study quantum computation in- 
volve manipulations of extremely mixed states. Recent 
research in error-correction and fault-tolerant computa- 
tion has shown that non-unitary evolution due to weak 
interactions with the environment results in no loss of 
computational power, if sufficiently pure states can be 
prepared Here we consider the situation where 

there are no errors or interactions with the environment, 
but the initial state is highly mixed. We investigate the 
power of one bit of quantum information available for 
computing, by which we mean that the input state is 
equivalent to having one bit in a pure state and arbitrar- 
ily many additional bits in a completely random state. 
The model of computation which consists of a classical 
computer with access to a state of this form is called de- 
terministic quantum computation with one quantum bit 
(DQC1). We demonstrate that in the presence of oracles, 
such a computer is less powerful than one with access to 
pure state bits. However, it can solve problems related to 
physics simulations for which no efficient classical algo- 
rithms are known. DQC1 is the first non-trivial entry in 
the class of models of computations which are between 
classical computation and standard quantum computa- 
tion. Investigations of such models are expected to lead 
to a better understanding of the reasons for the power of 



quantum computation. 

There are many kinds of problems that one might like 
to solve using a computational device. The three main 
problems not involving communication are function eval- 
uation, non-deterministic function evaluation and distri- 
bution sampling. Let S be the set of all binary strings 
and S n the set of binary strings of length n. The fc'th 
letter in a binary string b is denoted by bk- In function 
evaluation, we are given a function / : <S — ► S and a 
binary string b. The goal is to determine f(b). In non- 
deterministic function evaluation, we are given a relation 
R C S x iS and a binary string b. The goal is to determine 
a string c such that R(b, c). In distribution sampling, we 
are given a probability distribution P on S and a length 
n. The goal is return a string b of length n sampled 
according to P(b\n). 

In this paper we focus on deterministic function eval- 
uation. A one-bit function is a function whose values are 
single bits. Any deterministic function evaluation prob- 
lem reduces to that of evaluating a number of one-bit 
functions, one for each bit of the output string (a bound 
on the length of the output needs to be known). 

Any computation requires a certain amount of re- 
sources to obtain the desired answer. The most impor- 
tant resources are time and space. How much of either is 
required can depend on the computational power of the 
device, which determines the applicable model of compu- 
tation. For a given problem, one usually tries to deter- 
mine the resources required as a function of the problem 
size, which for one-bit functions is simply the number of 
bits in the input string. (For general functions, one can 
take it to be the sum of the sizes of the input and the 
output string.) The computational power of two models 
of computation is considered to be the same if for any 
algorithm in one model, there is a corresponding algo- 
rithm in the other one using at most a constant multiple 
of the resources of the first algorithm. An algorithm is 
considered to be efficient if the resource requirements 
are polynomial in the problem size. For the purpose of 
discussing efficient computation, the power of two mod- 
els is considered to be the same if for any algorithm in 
one model, there is an equivalent algorithm in the other 
model which uses at most a polynomial multiple of the 
resources. We adopt this definition for the remainder of 
the paper. In most cases below, the polynomial multiple 
is linear or constant. For a comprehensive treatment of 
classical computational complexity theory, see 0. 

The models of computation considered here subsume 
classical computation. The available computational de- 
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vices are assumed to include a classical probabilistic com- 
puter conforming (for example) to the model of an ab- 
stract random access machine with access to random bits, 
and a quantum system consisting of as many (quantum) 
bits as needed. Without loss of generality, we can assume 
that the procedure used to evaluate the function / on in- 
put b consists of repetitions of two parts. The first uses 
the classical computer to generate a sequence of unitary 
operations. The second applies these operations to the 
quantum system in a specified initial state. A measure- 
ment of the quantum system then yields the answer. An- 
swers from one repetition may be used to make decisions 
in the next ones. The models differ in what are the per- 
mitted initial states and measurements, and sometimes 
in the allowed unitary operations (a.k.a. quantum gates). 
In these models, function evaluation is performed prob- 
abilistically and we are satisfied with an answer which is 
correct with high probability, provided that we can make 
the probability of success high enough. 

The state space of the quantum system of n bits is 
the complex hilbert space Q n generated by the compu- 
tational basis. This basis is labeled by the strings in <S„; 
the basis element corresponding to string b is denoted by 
|6). Q n is the n-fold tensor product of Q 1 . To describe 
unitary operators and states, it is convenient to use the 
operator basis consisting of tensor products of the Pauli 
operators 
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The binary form of the indices consists of two-bit strings. 
A Pauli operator acting on the k'th bit is denoted by cr^ . 
A general tensor product of Pauli operators is denoted by 
Ob, where &2fc-i^2fc is the index of the operator acting on 
the fc'th bit. 

A pure state of the quantum system consisting of n 
bits is a unit vector in Q n . In general, the quantum 
system can be correlated with other systems that we do 
not have access to. As a result, the general state of the 
system can be described as a distribution of pure states. 
However, this description is overdetermined and a better 
one is based on the density operator formalism. A state 
is a positive definite trace one hermitian operator on Q™. 
Pure states are of the form (think of as a 

column vector, and ("01 as its hermitian transpose). All 
other states are obtained as convex combinations of pure 
states. 

When using highly mixed states, that is those with all 
eigenvalues small, it is convenient to describe the state by 



means of deviations from the identity expressed as sums 
of Pauli operators. In general, a state p can be written 
in the form 
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where the ab are real and is the bit string of all O's. 
The deviation of p is the non-identity component of the 
sum given by X^o afe(Tb - For example, the deviation of 
the state where the first bit is and the other bits are 
completely random is given by <Tz ■ 

The effect of an operation on the quantum system can 
be described by a unitary operator U which maps an 
input state to the output state U\ip). In terms of 
operators, U takes the state p to UpU^. 

In both of the models of quantum computation to be 
discussed in this letter, the elementary operations (quan- 
tum gates) that can be applied to the quantum system 
are a complete set of one and two bit unitary opera- 
tors ||,[l0) . Since global phase factors in operations have 
no effect on the computation, it suffices to have a com- 
plete set of determinant one operators. A simple such 
set consists of e" 10 ^^ for o^O (the tt/8 rotation) and 
the controlled tt rotation e ~^ k ?®^ 2 */4 g. This is a 
variant of the more familiar conditional sign-flip opera- 
tion. For fci = l,k 2 = 2, it takes |06) \Qb) 
and \lb) — ^"z 71 '/ 4 \\b). The difference in the effect on the 
second bit depending on the state of the first is a 7r ro- 
tation. More generally, one can use any U(b,t) = e~ l(7bt . 
The U(b, t) can be implemented with fidelity 1 — e using 
polylog(l/e)n of the previously mentioned set of oper- 
ations (this is a consequence of the results in jy]). A 
network implementation of a unitary operator U is a 
decomposition of U as a product of elementary gates, 

Deterministic quantum computation with pure 
states (DQCp): The initial state of the quantum sys- 
tem is |0) (where bold face means that all bits are in the 
state |0)). In the standard model of quantum compu- 
tation, the final (one-bit) answer is obtained by a mea- 
surement of the first bit. In DQCp, this measurement 
is replaced by a process which yields the noisy expec- 
tation of (Tz for the final state. For function evalua- 
tion (but not necessarily the other types of problems), 
there is no loss of power by making this restriction on 
the final measurement. To be specific, if the state of the 
quantum system is p, the measurement process returns a 
number which is sampled from a distribution with mean 
(ci^) = tr(cri 1 ^ p) and bounded variance s, where s is 
independent of the number of bits used. We call such a 
measurement an estimate of the quantity measured. Re- 
peating the computation and measurement are assumed 
to yield independent instances of this distribution. Thus, 
the mean can be estimated to within e with probability 
of error at most p using 0(log(l/p)/e 2 ) repetitions of 
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the computation |12[. DQCp is realized by an idealized 
bulk NMR quantum computer, where all molecules are 
initially perfectly polarized, there is no decoherence and 
there is no operational error. 

If U is a unitary operator with a network implementa- 
tion, then DQCp yields an estimate of tr(cri 1) £/|0)(0|[/ t ). 
By preprocessing the input state with operators of 
the form G = e~' lav7r ' 2 , any quantity of the form 
can be estimated. By further post- 
processing the state with one or two operators of the 
form G = e _4 °'<' lr / 4 , quantities like tr(cr c [/|b)(b|C/' 1 ') can 
be obtained. Since the resources required for pre- and 
postprocessing are linear in the number of bits, we can 
identify (the power of) DQCp with the ability to estimate 
these quantities. 

Deterministic quantum computation with one 
qubit (DQC1): The deviation of the initial state of the 
quantum system is The final answer is obtained as 
in DQCp by a bounded variance process yielding (o~ z )■ 
The initial state corresponds to having one bit in a pure 
state and the rest completely random. 

DQC1 is realized by an idealized high temperature 
NMR quantum computer where there is no decoherence 
and no operational error. This is the regime where the 
initial deviation state can be approximated by e~P H /Z ~ 

^r(I~PJ2i e i a ^) (noninteracting bits). The one-bit ini- 
tial state can be obtained by eliminating polarization in 
bits other than the first. Exploiting the additional po- 
larization can reduce the time resource required by at 
most a factor of 1/n, so that no additional power can 
be gained. Constraints on the amount of polarization 
cxtractable from the high temperature state have been 
investigated by Sorensen |L3]]. Note that DQC1 is not 
robust against many realistic error-models. For exam- 
ple, the results of Aharonov and Ben-Or jlj] show that 
error-correction cannot be used to successfully compen- 
sate against depolarizing independent errors due to the 
rapid loss of information if fresh ground-state bits cannot 
be introduced during the computation. 

Let U be a unitary operator with a network imple- 
mentation. Using DQC1 and pre- and postprocessing 
similar to that introduced with DQCp, we can estimate 
^{^aaUabW) for any a and b. We can therefore iden- 
tify DQC1 with the ability to perform these estimates. 

DQCp and DQC1 can be used to estimate any coef- 
ficient of certain operator expansions of U. For DQC1, 
write U — J2b ab<Tb ' with a b = wtr(<r;,Z7). This is the 
Pauli operator expansion of U. To determine a b , use the 
network for U to construct a network for the operator V 
which maps |0)|6) -> \0)U\b) and -> |1)|6). This 

is a "conditional" U operator and can be implemented 
with a linear amount of additional resources |l0| . Then 

^tr^Va^yt) = I(tr(C/ CTh ) + tr(a b U^)) (6) 



= Re(a b ), 

-^tr^Va^a^) = l {ttr (Ua b ) - itrfatft)) (7) 
= -Im(a b ). 

Thus the coefficients a b can be estimated with two com- 
putations. Since Ua b U' is a unitary operator easily im- 
plemented given networks for U, the ability to estimate 
the coefficients of the Pauli operator expansion is equal 
in power to DQC1. 

If the trick of the previous paragraph is used with 
DQCp, we can obtain any tr{U\b){b\). Since ti(U\a){b\) 
can be written as a sum of expressions of the form 
tr({7|(|o) + e^|6))((a| + e-^(b\), and the states in the 
second operator can be generated from |0) with 0(n) 
gates, any of the transition amplitudes of U can be deter- 
mined to get coefficients of the matrix U in the computa- 
tional basis. To see that the ability to estimate (a\U\b) is 
as powerful as DQCp, observe that tr(cri 1 ^[/|0)(0|J7 t ) = 
tr((0\U^aPu\0)). 

Let the evolution of a quantum mechanical system be 
described by a (possibly time varying) hamiltonian H on 
a hilbert space TL. To efficiently simulate this evolution 
using a quantum computer with n bits, it is sufficient 
to have a unitary embedding of TL in Q™ and an exten- 
sion H' of the embedded hamiltonian for which there 
are efficient quantum networks approximating e~ lH W s 
(taking h — 1) for small 6 to within 0(S 2 ). (To define 
efficiency properly requires problem instances for each n, 
in which case a resource bound polynomial in n and 1/8 
can be imposed.) A class of such hamiltonians consists 
of those which are a polynomially bounded sum of two 
bit hamiltonians j fl5| |. Whether it is useful to implement 
this simulation depends on what information one wants 
to obtain. If one simply wants to predict the outcome of 
a measurement in an experiment involving this system, it 
is possible to do that by simulation, provided that the ini- 
tial state can be computed, and the final measurement 
can be represented as a measurement in the computa- 
tional basis. DQCp and DQC1 are further restricted to 
measurements of expectations of operators which can be 
approximated as a sum of a reasonable number of com- 
putable conjugates of a z . 

Even if the initial state is restricted to the deviation 
cr^z \ there are no known classical algorithms for simulat- 
ing an arbitrary sum of two bit hamiltonians as described 
above. Since many real-world situations involve highly 
mixed initial states (e.g. most NMR experiments), such 
an algorithm is very interesting. In fact there are ongoing 
experiments exploiting the fact that many interactions 
can be simulated in solid state systems such as calcium- 
fluoride by modifying the dipolar interaction, but are dif- 
ficult to simulate using classical computation fl6[| . The 
observation that it it appears to be difficult to efficiently 
implement such simulations in DQC1 makes this model 
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useful and the question of the relationship between the 
powers of DQC1 and DQCp or classical computation 
non-trivial. 

For any hamiltonian, one of its physical properties of 
interest is the energy spectrum. This problem of deter- 
mining the spectral density function is not in the class 
of one-bit deterministic functions. However, any deci- 
sion problem based on the spectrum can be thought 
of as a non-deterministic one-bit function. The non- 
determinism is needed to permit the possibility that ei- 
ther decision is acceptable. Without this possibility there 
can be situations where we cannot gain arbitrarily high 
confidence in the correctness of the decision. 

In DQC1 it is possible to directly observe the spec- 
trum with a resolution inversely related to the effort used. 
Methods to accomplish this type of observation are well 
known in the NMR literature, where the natural hamil- 
tonians are routinely modified using effective hamilto- 
nian techniques fl7|j . Decoupling methods used to effec- 
tively remove interactions with some nuclei are an appli- 
cation |Q. Here we give a simple method for observing 
spectra from the point of view of quantum computation. 

Let U (t) — e~ lHt and assume that an efficient quantum 
network for applying U (t) with arbitrarily small error is 
available. Note that the quantum network may increase 
in complexity if less error is needed. Given the quantum 
network for U(t), we can derive networks for applying 
U(t) or U'(t) to bits 2 . . . n+ 1, conditionally on the state 
of the first bit. Let V(t) be the unitary operator which 
maps |1)|&) -> |l)tft(t/2)|&) and |0)|&) -» \0)U{t/2)\b). If 
we first apply a gate to transform the input state o* 
to cri , and apply V(t), then the deviation of the state 
becomes 

p, = ^(cos^yW + sin(A l i)4 1 ))| i )( i | (8) 

i 

= ^ J2 ( - C0S ( X ^ a - ) + sin(A;i)<T«) + rest. 

i 

where the Aj and |i) are a complete set of eigenvalues 
and corresponding eigenvectors (with repetition) of H. 
The second identity gives the expansion of pf in terms 
of (Tb's, with terms not of interest suppressed. The co- 
efficients of o x and o~ y can be measured and the re- 
sults combined into a single complex number with value 
f{t) = f^TT Si e ~ Ai *- This can be sampled at as many 
time points as desired and fourier transformed to obtain 
a broadened energy spectrum. The same technique can 
be used to measure the spectrum of a unitary operator by 
restricting t to integer multiples (corresponding to pow- 
ers of the operator). If the evolution V(t) is implemented 
as a power of V(At), and the measurement in DQC1 is 
non-destructive, then the spectrum can be observed di- 
rectly as is done in fourier transform NMR spectroscopy. 

DQC1 is strictly less powerful then DQCp in the pres- 
ence of oracles (a.k.a. black box). Suppose that we 



are given a black box implementing a unitary opera- 
tor U on our quantum system and we wish to deter- 
mine tr<0|t/ t CTi 1) C7|0), whose sign is the one bit answer 
computed by U on input |0). One method for obtain- 
ing this answer using DQC1 involves preparing a pseudo- 
pure state from the deviation Uz jl9[^0]|- Such a state 
has deviation proportional to that of a pure state. If the 
deviation is proportional to that of |0)(0|, measurement 

of after applying U is proportional to the desired 
answer. 

It turns out to be quite simple to prepare a kind of 
pseudo-pure state from using an ancilla bit, here 
labeled 0. Let T n be the unitary operator mapping 
|& 0) — > |(6 + 1)0) and for b^O, \b b) -► \b b) (addition 
of bits is modulo two). If we first swap bits and 1, then 
apply T n , and finally flip bit 0, the deviation is given by 
cri 0) (2|0)(0| -I). If we apply U to bits 1 through n, then 
the coefficient a of cr^cri 1 ^ in the state's deviation is the 
answer we want. As discussed above, this coefficient can 
be obtained using DQC1, with an intensity of a/2 n+1 . 
The exponential loss of intensity implies that it is very 
difficult to detect a above the noise. 

Unfortunately, without the ability to analyze a specifi- 
cation of U (for example, an implementation by a quan- 
tum network), we cannot do much better, even if we know 
that the answer of U is deterministic, that is a G { — 1,1}. 
The most general form of a DQC1 algorithm for deter- 
mining the answer can be described by k independent 
DQC1 computations consisting of quantum networks and 
calls to U and an inference function which attempts to 
determine a from the k measurements. Consider the first 
of these measurements. The expectation of the result of 
the measurement can be written in the form 

v(U) = ^tr(V r U...UV aMv3uK..UWJaU), (9) 

where r is the number of invocations of U, the Vi are the 
quantum networks used in the computation using m bits 
(of which m — n are ancillas). Since U is deterministic, 
f/|0)|O) = \b)\ip). If U is composed with the operator T 
which conditional on the state of bits 2 to n being 
flips the first bit, we get an operator U' which is also 
deterministic but has the opposite answer. We must be 
able to distinguish between the two. T can be written in 
the form T = I — 2P, where P is the pure state given by 
i(I— cr^)(g)|V')('0| (the m—n ancillas have been suppressed 
in this expression). For unitary operators W\ and Wi 
acting on m > n bits, we have 1 2^rtr(H / iPH / 2)| < 
By expanding U' = U — 2PU in the expression for v{U') 
we can write 

v(U>) (10) 

= ^r( tr (K(-2P){/y r _ 1 c/' . . .) + ti{v r uv r -iU' . . .)) 

= ^MV r (r-2P)UV r -iU'...) 
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+tr(V r UV r - 1 (-2P)UV r -2U' . . .) 
+tr(V r UVr-lUVr- 2 U' . . .)) 

^(tr(K(-2P)[/y r _it/'...) 

+tr(y r UV r - 1 (-2P)UV r - 2 U' ...) + ■•• + v(U)). 



Thus, w(£/') = ^-(ai 



«2r 



i>(t/)), where |ai| < 



2 m-n+i_ it follows that |u(lT')-u(l7)| < 4r / 2 "- It is clear 
that there is a sensitivity problem: To confidently dis- 
tinguish between U' and U requires exponentially many 
experiments or invocations of the oracle. 

The inequality derived in the previous paragraph sug- 
gests that it might be possible to obtain the answer of U 
with 0(2 n ) invocations of U rather than the 0(2 2 ™) ex- 
periments required to detect a signal of strength 1 /2™ in 
the presence of bounded noise. We do not know whether 
this is possible. 

We have introduced a new model of computation 
(DQC1) with power between classical computation and 
deterministic quantum computation with pure states 
(DQCp). DQC1 is not as powerful as DQCp in the pres- 
ence of unitary oracles but can solve problems in physics 
simulation for which there are no known efficient classi- 
cal algorithms. Since DQC1 requires only one quantum 
bit of information, the possibility that it adds power to 
classical computation is surprising. If this is the case, 
the usual reasons given for why quantum computation 
appears to be so powerful may have to be revised. On 
the other hand, a proof that DQC1 can be efficiently 
simulated by classical computation would be extremely 
interesting, as this could lead to practical algorithms for 
simulations of many experimentally interesting physical 
situations. 
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